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We consider a magnetic impurity in the antiferromag- 
netic spin-1/2 chain which is equivalent to the two-channel 
Kondo problem in terms of the field theoretical description. 
Using a modification of the transfer-matrix density matrix 
renormalization group (DMRG) we are able to determine lo- 
cal and global properties in the thermodynamic limit. The 
cross-over function for the impurity susceptibility is calcu- 
lated over a large temperature range, which exhibits univer- 
sal data-collapse. We are also able to determine the local 
susceptibilities near the impurity, which show an interesting 
competition of boundary effects. This results in quantitative 
predictions for experiments on doped spin-1/2 chains, which 
could observe two-channel Kondo physics directly. 



Magnetic impurities in low-dimensional antiferromag- 
nets are recently of much theoretical and experimental 
interest in connection with high temperature supercon- 
ductivity. We now study an impurity model in the spin- 
1/2 chain consisting of two altered bonds in the chain, 
which is known to have an equivalent field theory descrip- 
tion to the spin-sector of the two-channel Kondo (2CK) 
model (ly^]- The 2CK model has received much inter- 
est in the theoretical physics community 0] since it was 
first proposed in 1980 ||] and this model is often cited as 
a standard example of non-Fermi-liquid physics. Bethe 
ansatz || and conformal field theory techniques 0] have 
lead to almost complete understanding. On the other 
hand it is much less clear to what extend experimental 
applications exist, although there is some hope that the 
conductance behavior through metal constrictions can be 
explained by the 2CK effect Q . 

To study the 2CK effect experimentally, site-parity 
symmetric bond defects could feasibly be created by 
doping quasi one-dimensional spin-1/2 compounds, and 
Knight shift experiments would then be able to observe 
2CK physics explicitly. Our goal is therefore to provide 
quantitative predictions for the local susceptibilities in 
a range around the impurity, which turn out to exhibit 
an interesting competition between two boundary effects. 
Moreover, we are able to explicitly show the expected 
data collapse of the impurity susceptibility for different 
coupling strengths, and calculate the resulting cross-over 
function over a large temperature range. For this purpose 
we have developed a modification of the transfer-matrix 
DMRG @. 

The model we are considering is the antiferromagnetic 



spin-1/2 chain with two altered bonds 



L-l 



(1) 



which is known to have an equivalent field theory descrip- 
tion of the 2CK effect as we will outline below. Interest- 
ingly, an integrable spin-1 chain with a spin-1/2 impurity 
is also known to be equivalent to the 2CK problem on the 
level of the Bethe ansatz equations [||, but we arc not 
able to analyze that model in terms of the field theory. 

The low-temperature and long-wavelength properties 
of the unperturbed spin-1/2 chain are well described by 
a conformal field theory Hamiltonian in terms of 1+1 
dimensional bosons 



H 



dx - [{d^f 



(2) 



where is the conjugate operator to the boson field <f>. 
This field theory has been discussed in more detail else- 
where and we will just focus on a more pedagogical 
description in this Letter. For temperatures and wave- 
vectors well below some cutoff A ~ J the spin-1/2 chain 
exhibits critical behavior and scale invariance. Scale- 
invariance means that we get the same physical results 
for some quantity O if we rescale the temperature or dis- 
tances with some factor T, as long as we also multiply 
the physical quantity O with r d , 



0{T) = T- d O{TT), 



(3) 



where d is referred to as the scaling dimension of 0, and 
the temperature T can also be replaced by the inverse 
system size 1/L. For example, operators in the "free" 
Hamiltonian (^) have a scaling dimension of d = 1 (after 
integration over x). This means that the energy spacing 
of the spectrum is proportional to 1/L. "Higher order 
operators" in the Hamiltonian are operators with d > 1, 
which give small corrections to the spectrum of order 
1 / L d and are hence termed "irrelevant" and are neglected 
in Eq 1 ph . (We have also neglected one marginal operator 
cos ^/8tt(/) which turns out to be justified although the 
corrections are only logarithmically small). For the rest 
of this Letter we will work in the thermodynamic limit 
L — > oo and rescale T instead. 

We can now analyze the perturbation J' on the bonds 
in Eq. (Q) in terms of the scaling dimensions of the lo- 
cal operators, which arise in Eq. (j^) due to the broken 
translational invariance. If the perturbation on the bonds 
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5 J = J— J' is small, we can use perturbation theory on a 
chain with periodic boundary conditions. In this case the 
local operator with the lowest scaling dimension which 
still observes the site-parity symmetry of the problem is 
known to be d x sin of dimension d = 3/2 (lj. This 

operator gives the leading corrections, but is still irrele- 
vant. In particular, for a given coupling strength S J the 
size of the corrections becomes effectively smaller by a 
factor of r d_1 if we rescale T by T < 1. The opposite 
is also true: If we increase the temperature the effec- 
tive perturbation strength r d-1 <5 J may become so strong 
that a systematic expansion fails at some special temper- 
ature Tk, called the cross-over (or Kondo) temperature. 
Above Tk we therefore expect a completely different be- 
havior, namely that of an open chain if Tk < T < A. We 
say that the system renormalizes from an open bound- 
ary condition to the infrared fixed point of a healed pe- 
riodic chain as the temperature is lowered. For small 
5 J we have defined Tk as the temperature at which the 
product T^SJ becomes large, so that we can write 
Tk oc (5J 1 ' = 6J~ 2 in this limit. Hence, rescaling 
T by r is equivalent to changing Tk by 1/r, i.e. altering 
the initial coupling strength 5 J by r^ 1 = T 1 / 2 , which is 
really the meaning of renormalization. 

To compare this system with the 2CK effect it is more 
instructive to start with open boundary conditions and 
consider a weak antiferromagnetic coupling < J' -C J. 
In this case the leading operator is S§[d x cf)(Q) + d x <f>(L)] 
with dimension d = 1 which turns out to give logarith- 
mically relevant contributions as the temperature is low- 
ered. For small J', the Kondo-temperature is given by 
Tk oc e~ b / J , where b is some constant. If we identify 
the central spin Sq with the Kondo-impurity and the two 
ends of the chain with the spin-sectors of the two elec- 
tron channels, it is evident how this scenario is equivalent 
to the 2CK problem f|: A small coupling (J' <C J) to 
the impurity is marginally relevant as we lower the tem- 
perature and the system renormalizes to an intermediate 
coupling (J' = J) fixed point, i.e. the periodic chain. 

We have chosen the numerical transfer matrix DMRG 
jjj to test these concepts. The partition function of the 
spin-1/2 chain can be written in terms of transfer ma- 
trices Z = tr(T L / 2 ) A L / 2 , where A is the largest 
eigenvalue of T. The transfer matrix DMRG constructs 
a transfer matrix for a finite number of time-slices M and 
then successively increases M by keeping only the most 
relevant basis states that are necessary to calculate A. 
We now extend this method to non-uniform systems by 
using the eigenstate \ip\) corresponding to the highest 
eigenvalue A. We can include any impurity interaction 
which is described by a local matrix 7i oca i explicitly in 
the partition function 

Z = tx{T L / 2 ~ l T loca {) A^-^Al^oeall^). (4) 

If the eigenstate \ip\) is known to reliable precision from 



the DMRG, it is straightforward to determine any ther- 
modynamic property from this expression (even locally) . 
This method proved to be superior in speed, accuracy, 
and temperature range compared to Monte Carlo simu- 
lations, which we have used for checking purposes. 
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FIG. 1. The linear response of the impurity spin to a lo- 
cal magnetic field. Inset: the data-collapse for 14 different 
J' = O.IJ, J and the effect of a finite field Bo = A c h. 

The first task is to show that our renormalization pic- 
ture is accurate i.e. that we can indeed use periodic 
boundary conditions to describe the system at sufficiently 
low temperatures. For this purpose we consider the linear 
response xo of the impurity spin Sq to a local magnetic 
field Bq, which is given by the Kubo formula 

X0 (T) = (S^(r)S^O)) dr. (5) 

For open boundary conditions the impurity spin will 
be asymptotically free and the response is given by a 
Curie law xo oc I/4T. On the other hand, for peri- 
odic boundary conditions the auto-correlation functions 
will be determined by the leading field theory operator 
which obeys the same symmetry transformations as Sq. 
We found this operator to be cos-s/^t/) with dimension 
d = 1/2 and an auto-correlation function proportional 
to f/r, which should lead to a logarithmic divergence 
as T — > 0. As can be seen in Fig. (|f|), we indeed find 
a cross-over from Curie law to logarithmic behavior at 
an effective cutoff A e g = min(A, Tk), depending on the 
coupling strength J' with a scaling behavior of the form 
XoCO = g(T/A e f[)/A e ff (inset). The logarithmic behav- 
ior was observed before at one coupling only by other 
methods Note, however, that the response xo is 

only an indication of the auto-correlation functions, but 
it must not be confused with the impurity susceptibil- 
ity, since we have only applied the magnetic field at one 



2 



spin. In fact every spin in the periodic chain shows a 
logarithmic response to a local magnetic field. 
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FIG. 2. The scaled impurity susceptibility TxXimp for an 



appropriate choice of Tk as a function of J' 
(inset). 
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From an experimental point of view it is much more 
interesting to look at the impurity susceptibility, which 
can be defined as the size independent contribution to the 



total system susceptibility \ 



imp 



linir 



L — >oo (Xsystcm" 



■Lx), 

where x is the susceptibility per site far away from any 
boundary. While our method is in principle capable of 
extracting the impurity susceptibility directly, we were 
able to obtain results at much lower temperatures by 
considering the linear response of the impurity spin to a 
uniform field B instead 



Xh 



d(SS)B 
dB 



J'x, 



(6) 



which gives a good indication of the true impurity sus- 
ceptibility (in contrast to xo in Eq. (|J), where we only 
considered a local field Bo). From the 2CK effect it is 
known that the impurity susceptibility is logarithmically 
divergent below Tk while it shows a Curie-law behav- 
ior above ||||- Interestingly, this cross-over shows a 
universal data collapse, because changing the coupling 
strength (i.e. Tk) is equivalent to rescaling the tempera- 
ture, i.e. there is only one independent variable T/Tk 



WT) = f(T/T K )/Ti 



K ■ 



(7) 



where f{x) is a universal function (ignoring higher order 
operators). This data collapse is clearly seen in Fig. (^) 
with an appropriate choice of Tk as a function of J' (in- 
set), showing the predicted logarithmic scaling at low 
T . The non-universal deviations of some of the curves 
at higher T/Tk are due to the fact that in those cases 



Tk was so large that regions above the cutoff have been 
included. 

In principle, a similar logarithmic scaling should be ob- 
servable for the impurity specific heat C lmp /T. However, 
we find that the critical scaling of the specific heat occurs 
at lower T and instead Ci mp shows a more complex be- 
havior in the intermediate range. Moreover, the numeri- 
cal method is known to produce larger errors for C at low 
T , which may be due to the second derivative involved. 
Therefore, we could not fully reproduce the correspond- 
ing data collapse of C- imp /T, but our data is nonethe- 
less consistent with a cross-over to logarithmic scaling as 
T — > 0. Just like in the 2CK effect we also find that a 
finite magnetic field (local or global) will change the log- 
arithmic non-Fermi-liquid behavior and produce a cross- 
over to a constant susceptibility below some temperature 
Tfl as shown in inset of Fig. (|l|). The different renor- 
malization behavior can be traced to the broken spin-flip 
symmetry, which allows the relevant operator cos \/2tt<P 
with d = 1/2 at the periodic chain fixed point. This will 
result in a quadratic field dependence Tfl oc Bq, which 
has already been demonstrated in Ref. B. In experi- 
ments, however, fields are expected to be small compared 
to J and Tk, so that in most cases only the non- Fermi 
liquid behavior will be observed. The effect of a finite Bq 
is analogous to violating the symmetry condition in the 
resonant tunneling scenario for electrons pj . 

Finally, we would like to consider the local susceptibil- 
ities (the Knight shift) of the individual spin sites in a 
region around the impurity as a function of site index x 



XlocalO) 



d{S z {x)) l 



dB 



(8) 



B=0 



In Ref. [[Tl| the most dramatic effect of an open boundary 
condition on the Knight shifts was found to be a staggered 
component Xopcn = Xiocai^ X m the response to a uniform 
magnetic field B, which actually increases with site index 



Xopen(lE) OC (-1) 3 



xVt 



V sinh AxT 



(9) 



where T is measured in units of J. 

The staggered part in Eq. @ arises due to open bound- 
ary conditions and hence it will be diminished as the sys- 
tem renormalizes to the periodic chain fixed point. How- 
ever, there will be a whole new effect due to the lead- 
ing irrelevant operator d x sin y/2n^> at the periodic chain 
fixed point. This operator also induces a staggered part 
Xpcriodic, but with opposite sign, i.e. the induced response 
at the first site Si is negative. The alternating response 
as a function of site index x is now 

Xpcriodicfz) = (-l) 2 ^ Jdy {S* lt (x)S* ni (y)) 

cxf-lfi fdy f 1 drg(x,y,r), (10) 
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where S^L and S* ni refer to the leading operators which 



describe the alternating (cosv2tt0) and uniform (d x 4>) 
parts of the spin z-component, respectively. The correla- 
tion function g(x,y,r) is therefore given by 



g(x,y,r) oc (cosV27rc6(x) d^y) d x sinV 27rc/>(0, r)) 
T T 2 



OC 



sinT(7TT - i2x) sin 2 T(itt - i2y) 



(11) 



where we have used standard field theory techniques 
(e.g. the boson mode expansion E3]). The integral of the 
second factor over the spatial coordinate y is the same 
that determines the unperturbed susceptibility per site 
|l3f and simply gives a r-independent contribution pro- 
portional to T. The integral of the first factor over the 
imaginary time variable can then easily be done, giving 



Xpcriodic(^) oc (-1)' T log [tanh(xT)]. 



(12) 



This expression shows the logarithmic divergence with T 
explicitly for small x at the impurity, and the response 
then drops off with exp (— 2xT) as xT — > oo (i.e. with the 
same exponential as observed for the open chain response 
(§)■ 

As J' is increased, the alternating part changes from 
the behavior of Eq. (||) to the behavior of the stable fixed 
point in Eq. (|l2]), which is always logarithmically domi- 
nant as T — > 0. However, even very close to the periodic 
chain fixed point we observe an interesting and complex 
competition of both contributions (see Fig. (^a)), which 
is nonetheless completely understood. In particular, be- 
low Tk the total amplitude of the staggered part of xiocai 
always fits very well to a superposition 



Xtotex(x) = ci log[tanh(xT)] + c 2 



xVT 



V sinh 4xT 



(13) 



This formula gave excellent results as can be seen for 
a typical fit in Fig. (|]c), and the coefficients have been 
determined for all values of J' > 0.2 J and T (inset). The 
coefficient c± is T-independent, while c-i renormalizes to 
zero as T — ► 0. 

To see these effects experimentally, it will be necessary 
to induce site-symmetric bond defects into quasi one- 
dimensional spin- 1/2 compounds like Sr2Cu03, KCUF3, 
or Copper Pyrazine Nitrate. A simple doping with an- 
other effective spin-1/2 ion for Cu 2+ may be possible, 
but more likely the surrounding non-magnetic ions can 
be doped in order to create a suitable lattice deformation. 
The effect of open boundaries has already been seen in 
NMR experiments [H, by observing a unique feature 
that broadened with 1 / y/T and had the predicted shape 
derived from Eq. (H). For site-symmetric perturbations 
we can now predict a broadening with log T of the NMR 
spectrum. In addition there will be a feature (kink) at 
smaller Knight shifts which comes from the relative max- 
imum in the alternating response (see Fig. (|b)). This 



feature has a distinctive shape and temperature depen- 
dence which can be calculated from the constants c\ and 
C2 with the help of Eq. (|l3|) for any coupling strength 
J'. Ordinary susceptibility measurements should be able 
to identify the impurity contribution as predicted by the 
cross-over function in Fig. (Q). In a more exotic twist 
we can even imagine /iSR experiments, where the muon 
itself may play the role of an impurity, i.e. depending 
on the preferred muon location in the lattice, the muon 
could feasibly induce a site-symmetric lattice distortion 
and at the same time measure the logarithmically diver- 
gent response at the impurity site. 

In summary, we have shown the cross-over of a two- 
channel Kondo impurity explicitly, which confirmed the 
renormalization picture and demonstrated the expected 
data collapse explicitly. Moreover, we were able to pre- 
dict the response to a uniform magnetic field of individual 
spin sites in a large region around the impurity, which led 
to quantitative predictions for Knight shift experiments 
on doped spin-1/2 compounds. 
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FIG. 3. a.) The local susceptibilities as a function of site 
index for T = 0.05J and J' = 0.7J. b) The corresponding 
typical NMR spectrum with a distinct feature (kink), c) The 
fit of the alternating amplitude to Eq. ([L3|) (solid line) with 
the appropriate coefficients (inset). 
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